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In this paper we analyze the ground state of the heavy-quark qqqG system using standard princi- 
ples of quark confinement and massive constituent gluons as established in the center-vortex picture. 
The known string tension Kf and approximately-known gluon mass M lead to a precise specifica- 
tion of the long-range non-relativistic part of the potential binding the gluon to the quarks with 
no undetermined phenomenological parameters, in the limit of large interquark separation R. Our 
major tool (also used earlier by Simonov) is the use of proper-time methods to describe gluon prop- 
agation within the quark system, along with some elementary group theory describing the gluon 
Wilson line as a composite of co-located q and q lines. We show that (aside from color-Coulomb 
and similar terms) the gluon potential energy in the presence of quarks is accurately described (for 
small gluon fluctuations) via attaching these three strings to the gluon, which in equilibrium sits 
at the Steiner point of the Y-shaped string network joining the three quarks. The gluon under- 
goes small harmonic fluctuations that slightly stretch these strings and quasi-confine the gluon to 
the neighborhood of the Steiner point. To describe non-relativistic ground-state gluonic fiuctua- 
tions at large R we use the Schrodinger equation, ignoring mixing with / = 2 states. Available 
lattice data and real-world hybrids require consideration of R values small enough for significant 
relativistic corrections, which we apply using a variational principle for the relativistic harmonic 
oscillator. We also consider the role of color-Coulomb contributions. In terms of interquark sepa- 
rations R, we find leading non-relativistic large- J? terms in the gluon excitation energy of the form 
e{R) -> M + £,[Kf/{MR)Y^'^ - (ac/R where are calculable numerical coefficients and Qc — 
0.15 is the color-Coulomb qq coupling. When the gluon is relativistic, e ~ {Kf / R)^^^ ■ We get an 
acceptable fit to lattice data with M = 500 MeV. Although we do not consider it in full detail, we 
show that in the qqG hybrid the gluon is a bead that can slide without friction on a string joining 
the q and q. We comment briefly on the significance of our findings to fluctuations of the minimal 
surface, a subject difficult to understand from the point of view of center vortices. 

PACS numbers: 11.15.-q, 12.38.-t, ll.15.Tk UCLA/04/TEP/51 



I. INTRODUCTION 



In QCD, hybrids arc hadrons that cannot be understood as combinations of quarks and antiquarks only, but that 
require consideration of gluonic excitations along with the quarks. Perhaps the most interesting of these are the exotic 
hybrids, which are hadrons whose quantum numbers are not found in any colorless combination of q and q, and so 
must contain valence gluons. Hybrids and exotics jj have been studied theoretically for decades (for example, 0], 
(3), but these early studies are tied to models and approximations that hardly make it possible to choose one model 
from another. Experiments have searched for hybrids as well, and there are a number of credible observed candidates 
for mesonic hybrids (for a recent review see The author knows of no plausible candidates for baryonic hybrids. 

Unfortunately, theoretical approaches to hybrids are, naturally enough, even more forbidding than approaches to 
ordinary hadrons. So one might wonder what the value of a theoretical paper on hybrids is. The answer is that 
recent lattice simulations 0, of heavy static quarks (described by straight Wilson lines in the Euclidean time 
direction) with single-gluon excitations provide a valuable testing-ground for theoretical ideas, because one does not 
have to deal with the extra complication of quark Wilson-line fluctuations and because lattice data are available at 
theoretically- interesting separations of quarks. 

We will be primarily interested in estimating the ground-state gluonic excitation energy e{R) of the heavy-quark 
baryonic hybrid qqqG, as a function of a suitably-defined interquark separation R, and in comparing our results to 
lattice data 0. We will also remark in less detail on similar issues for mesonic hybrids, but will not attempt to 
compare to data in this paper. The regime where we are on firmest ground is the regime of asymptotically-large large 
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R, where the dynamics is non-relativistic (because the gluon has a dynamical mass) and we can fairly confidently 
use the Schrodinger equation. However, lattice data are available for smaller R, where relativistic corrections are 
important, and to cover the entire relevant range of R we use a variational principle incorporating fuUy-relativistic 
kinetic energy for the giuon. The smallness parameter for non-relativistic motion is a numerical constant times 
[Kf/{M'^R)Y^^ , where Kp is the string tension and R the intcrquark separation. It turns out that our calculations 
are in excellent agreement with the lattice data for a gluon mass of 500 MeV, even for R small enough for relativity 
to matter. 

One should note that in the fuUy-relativistic regime there should be mixing with other states containing more 
gluons, such as baryon plus glueball. Additionally, the derived potentials have angular momentum I = 2 terms as well 
as ^ = 0. These lead to mixing of gluonic I = 2 and I = states, and should lower the energy slightly compared to the 
results we give. This lowering might be compensated by raising the gluon mass. We do not consider either of these 
mixings here. 

The idea that the gluon has a constituent mass M 0| goes back a long way, but a purely theoretical determination 
of M is not yet within our powers. Lattice simulations have verified this idea and found values for the mass, although 
many such simulations, done in particular gauges, cannot guarantee a gauge-invariant determination of the mass. 
However, Ref. 8] has extrapolated their Euclidean lattice data to the Minkowski regime, to find the so-called pole 
mass, which is gauge-invariant. The pole mass so found on the lattice is M ~ 600 (-f 150 -30) MeV We believe 
that our fitted value of 500 MeV is acceptably close to the lattice pole mass, given the expected uncertainties in our 
approach and in the lattice determinations. 

There are two main steps in finding the dependence of the gluonic energy on R. The first uses a proper-time 
technique for describing propagation of the gluon in the hybrid, and the second, as described above, implements the 
dynamical fluctuations in the proper-time gluon paths with a Schrodinger equation modified for relativistic effects. 
In the first step, we describe the propagator of the gluon in the heavy-quark hybrid as a quantum perturbation on 
a background of confining gauge potentials. The center-vortex picture of confinenment prescribes llOj l that the 
quarks are joined by a Y-shaped network of minimal surfaces (or strings, at fixed time). The gluon propagator has a 
proper-time form that is an integral over all paths of an adjoint Wilson line going along a path joined at both ends 
to the system of heavy-quark Wilson lines. [These techniques were developed earlier by Simonov [Tlj and applied 
to qqG hybrids, but not qqqG hybrids. Simonov's works were unknown to the present author, who developed them 
independently and discovered Simonov's work after this paper was finished.] As we will see, the result is that if we do 
not distinguish the three quarks by flavor or other quantum numbers, the gluon can be described as a massive bead 
attached to the three strings joining the quarks in the hybrid, and fluctuating in directions transverse to the string. 
The fluctuations stretch the string slightly, and there is a harmonic restoring force on the gluonic bead. In the qqG 
hybrid the gluon is a bead sliding without friction on the string joining the quark and antiquark. 

All the potentials we use in the second step are harmonic, found by expansion of various minimal areas occurring in 
the baryonic Wilson loop. Except in the ultra-relativistic limit and the Newtonian limit this relativistic Schrodinger 
equation is not analytically solvable, but good approximations to s{R) for large R can be found by a variational 
technique. We describe this technique in the Appendix. For the qqqG system, the color-Coulomb energy can be 
straightforwardly estimated classically and added as a perturbation. It is not quite so straightforward for the qqG 
hybrid color-Coulomb energy, which we describe in the Appendix. The Appendix also treats the string energy for the 
qqG hybrid, flnding relativistic gluon energies whose leading term is of the form Nn/ R. Lattice data on heavy-quark 
qqG systems suggest such leading string-like energies, but with substantial corrections 6, 12]. These corrections can 
be estimated with our techniques, but we do not study the problem here. 

Our final result yields three terms in the large- i? (non-relativistic) expansion of the qqqG energy: 

where ^, ^ are numerical parameters that we estimate with fair accuracy below and in the Appendix, and ac is the 
coefficient of — in the qq Coulomb energy. At shorter distances, where relativistic corrections are important, 
the gluon energy scales like {Kp/RY^^ , and Coulomb corrections are smeared out by the gluon wave function, as 
discussed in the Appendix. 

There are numerous other corrections to the small- i? results, including spin-dependent terms and poorly-understood 
terms that arise when two Wilson loops partly coincide, as is the case for quark Wilson loops formed from the 
gluon propagator. We ignore all such terms, which do not seem to be very large, judged by our fit to the lattice 
data. It is somewhat surprising that our techniques give good results even for fairly short distances, considering all 
the possible effects that can enter. 
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II. PROPER-TIME TECHNIQUES AND WILSON LOOPS 

In gauge-theory dynamics, a Wilson line or loop is to be integrated not only over all gauge configurations but also 
over all possible lines joining the fixed endpoints, with a certain weight function. This weight function encodes all 
the particulars about the mass, flavor, and spin of the particle being described by the Wilson line. In other words, 
particle propagation is to be described by a proper-time path integral. 

A. Gluon proper-time propagator 

We wish to construct the necessary path integral for a gluon in a hybrid state, where distance scales are large 
compared to the QCD scale. We will argue below that gluonic spin contributes only short-range effects, which allows 
us to simplify matters by considering the propagation of a scalar particle in the adjoint representation. In a fixed 
background gauge potential the Euclidean proper-time propagator A(a;; y) from x to y oi a scalar particle of mass M 
is 

Aix;y)=^ / (dz)exp{— / dT[P + l]}UA{x;y) (2) 

Jo Jx ^ Jo 

where (dz) stands for the integral over all paths ^^(r) from x at t~0 to y at t = s; A/" is a normalization factor; and 
Ua{x; y) is the adjoint Wilson line 

UA{x;y) = PeMw J dTi^T"A^(z)]. (3) 

Here T" are the adjoint generators, and P stands for path-ordering. Of course, the propagator A is not gauge-invariant 
by itself, but we will soon couple it to quarks in a gauge-invariant way. 

If (still ignoring spin) we identify A as the propagator of the potential A'^{z) that occurs in Ua, intractable problems 
arise. So we envision the gauge potential as being separated into a confining part (condensate of center vortices) that 
constitutes the background potential and the small-amplitude part whose propagation we are describing; it is this 
latter part that constitutes the gluonic excitation in a hybrid hadron. As for spin itself, the propagation of a vector 
gauge boson with small amplitude in a background field requires that in Ua be replaced by Zf^ + Y,^iyd^, where Y^j^^dn 
are the spin-one generators of rotations in Euclidean four-space, and path ordering is also carried out with respect 
to the spin operator. However, the derivative yields the background potential field strength, which is short-ranged. 
Since our primary interest is in long-range effects, we will from now on ignore gluon spin, and think in terms of the 
explicit scalar propagator of Eq. Q . 

B. Splitting the adjoint Wilson line 

It is familiar fact that a gluon line in a Feynman diagram can be replaced by co-located quark and antiquark lines, 
joined to other quarks according to certain rules. The usual large- A'^ counting rules, for example, follow from such 
considerations. Perhaps less familiar is the dynamical application of this splitting of a gluon line ,,1L |. which we now 
describe. The resulting picture is most easily described for qqG hybrids, which we discuss first. 

1. Mesonic hybrids 

Let 

H,{x) = V^74^"V^G^-(^) (4) 

be the interpolating field for a hybrid meson, with (l/2)Act the conventionally-normalized generator in the fundamental 
representation. (It happens to be the field describing an exotic, with J^^ — 1 Q.) The field ip describes very 
massive quarks. Ignoring spin indices and other irrelevant complications, the H-6.e\d propagator Ah is an expectation 
value 



AH^{Uiil)U^{2)[UA{9)]li) 



(5) 




FIG. 1: Left: A rectangular Wilson loop for a quark (1) and antiquark (2) and a gluon (G). Right: The same loop with the 
gluon line decomposed into another quark and antiquark line, forming two standard Wilson loops. 



where 17(1,2) are fundamental Wilson lines for the quark and antiquark, and (•) means to integrate over all gauge 
potentials and over the gluon Wilson line as shown in Eq. (|2Jl . The gluon Wilson line matrix has been converted into 
a matrix with fundamental indices on it, via 

[UAi9)]ii^\i\a)U^p)i[UAi9)U. (6) 

Here Ua is the adjoint representative of the group element g, which is the path-ordered product seen in Eq. ||2Jl. We 
show the index structure of Eqs. in Fig. ^ 

Now we split the gluon line, using the SU{N) equations 

li>^a)U^p)i^SiSl-^Si6i (7) 

and 

[Ua{9)U = lTr[U{g)X^UHg)Xp]. (8) 

Combining leads to a familiar equation: 

[UA{9)t = Ul{g)[U\g)]i-^SfSi. (9) 
This equation has the group property: If <? = 5i.92, the product 

[UA{gi)]inUA{92)ti, (10) 

when computed using the right-hand side of Eq. @, has the form of that right-hand side with g replaced by gig2- 

For us, g is the group element in the gluon Wilson line [see Eq. ^]. The adjoint on the right-hand side of Eq. © 
means to trace the path backwards, as one easily checks. Therefore the hybrid propagator in Eq. Q assumes the 
simple form 

Ah{x; y) = {TrU{lG)TrU{G2)) + ... (11) 

where IG is the closed fundamental Wilson loop shown in the right side of Fig. ^ and similarly for the other closed 
loop G2. The omitted term comes from the 1/7V term in Eq. lO, and yields a single Wilson loop describing a qq pair. 
We will drop this term, which in any case cannot contribute for exotic hybrids. 
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FIG. 2: Left: Wilson loop for a qqqG configuration (quarks labeled 1, 2, 3; gluon labeled G). Right: The same loop with the 
gluon line decomposed into another quark and antiquark line, with quark line 3 singled out as discussed in the text. 

2. Baryonic hybrids 

Similar considerations apply for qqqG hybrids, but there are minor group-theoretic complications. The direct 
product of three quark representations is 

3(^3(g)3 = le8®8'el0 (12) 

There are two 8s which can couple to a gluon. Denote the quark group vectors by a% 6*^,0', and form the octet 
combination 

V;{a- b, c) = e.kia'h'^c' - i(5jep,,aP6«c^ (13) 

Quark a has been singled out (its index is external). By permuting the quark group vectors one seems to come to 
three possible octets, but because of 

V;{a- b, c) + VJib; a, c) + Vjic; a,b) = (14) 

there are really only two. 

We can construct a propagator for a qqqG hybrid in the same spirit as the earlier work for the mesonic hybrid. Let 
there be three heavy-quark lines labeled 1, 2, and 3, plus gluon line G. We write the propagator for one of the two 
octets, singling out line 3 as we did quark a above (again dropping the trace term, which is missing the gluon): 

Using the identity of Eq. @ as before yields the qqqG analog of the gluon line splitting for mesons as shown graphically 
in Fig. 13 Line 3, which has been singled out, forms a qq state with one of the split gluon lines, while the other split 
gluon line forms a new baryon 12G. 

C. Gluon dynamics 

In the quark Wilson loops of Figs. the rectangular lines do not fluctuate (the quarks are infinitely heavy). But 
the lines associated with the gluon do fluctuate, since the gluon path-integral weight function has a finite mass. The 
dynamics of the quark loops is governed by conventional area laws, and we now turn to the question of how those 
area laws govern the gluon dynamics. 



FIG. 3: An adjoint Wilson loop displayed as the superposition of two oppositely-oriented Wilson loops. 



The center- vortex picture of the area laws for a fundamental Wilson loop claims |2J that, without the extra gluonic 
effects of the present paper, the area to be used in the area law is the minimal area spanning the given Wilson loop. 
(This only makes sense if there is a unique such area, which is certainly true in the present case; see 9] for further 
discussion.) 

D. Qualitative behavior of gluonic potentials 

Let us begin with some qualitative remarks. Consider a Wilson loop in the adjoint representation, which we 
decompose into two fundamental (quark) Wilson loops in the now-familiar way as shown in Fig. |31 Because the 
superposed loops are oppositely-oriented, the linkage of a center vortex with one loop is exactly and coherently 
cancelled by the linkage of the same vortex to the other (oppositely-oriented) loop. Thus the center-vortex linkage 
contributes a factor of unity to the adjoint Wilson loop, as must be the case. 

Consider the right-hand part of Fig. ^ showing the gluon world line decomposed as a quark and antiquark line. 
The area law for the original qq pair, labeled 1 and 2, is the sum of the two areas on the right (which of course is 
the original area on the left side of the figure), as long as the gluon line stays in the original minimal surface (plane 
of the figure) and does not stray too close to the lines 1 and 2. (As the gluon line approaches these lines, perimeter 
terms arise because one or the other of the two Wilson loops in the right-hand part of the figure has a length scale 
comparable to the QCD scale). That is, fluctuations of the gluon line which are confined to the original minimal 
surface change nothing; the gluon does not feel the qq pair. 

If the gluon line fluctuates to a position not on the original minimal surface, there are forces acting to pull it back. 
Just as for the adjoint Wilson loop of Fig. O these are not confining potentials (area laws), but are due to perimeter 
and other forces, including color-Coulomb forces. Consider first a gross deformation of the Wilson loop of Fig. 
(right), as shown in Fig. 0] Although the areas enclosed by contours 1 and 2 have increased, just as for Fig. ^they 
are oppositely-oriented. The only area law refers to the original 1-2 loop. In fact, if one tries to pull the gluon a long 
way from the qq loop, the adjoint string discussed above will break, forming a gluon pair, and the result will be a 
hybrid plus a glueball. 

Evidently the techniques described here are applicable to the problem of fluctuations of the area in a confining area 
law, although this fluctuation problem is more complicated because one must treat an arbitrary number of gluons, 
not part of the condensate of center vortices, exchanged between parts of the Wilson line. Originally, area-law surface 
fluctuations were thought of as taking place in chromoelectric flux tubes ;13.] . However, the center vortex description 
of confinement is different, and requires a new approach. We discuss this a little more in the concluding section. 
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FIG. 4: The gluon, decomposed as a qq pair, is pulled far from the original minimal surface of Fig. Q 



E. Toward quantitative gluonic potentials 



If the two Wilson loops IG and G2 shown in the right side of Fig. ^were statistically independent with respect 
to the center vortices linking them, the expectation value in Eq. p II for the hybrid propagator would factor, and we 
would have for this propagator 



^h{x-v) ^ N / (dz)exp{— / dT[P + l]}{TrU{lG))A{TrU{G2))A 

JQ Jx ^ Jo 

poo py ps 

= Af / (dz)exp{—- dT[P + I]}exp[--KF{A{1G)+A{G2)] 
Jo Jx 2 Jq 



(16) 



where {■)a is the expectation value over center vortex configurations and we have explicitly exhibited the gluonic path 
integral. In this equation A(1G) is the minimal surface spanning contour IG. Actually, the center vortex linkages in 
the two loops are not quite independent, because the loops share the leg G, but the correlation effects are short-ranged 
(see |j for a discussion) and will be ignored here. 

Take the original Wilson loop 12 of Fig. ^ to be a rectangle of length T in the Euclidean time direction t, and of 
length R in the z direction, centered at 2 = 0. We will extract the non-relativistic potential from the path integral of 
Eq. (|16|l as the coefficient of the length T of the time-directed segments of the original quarks paths 1,2, in the limit 
— >■ 00. As mentioned previously, all potentials for hybrids will be expressed relative to the potential in the original qq 
system, so we subtract RT from the area ^(IG) -I- A{G2) in determining the gluon potential. If the gluonic contour 
segment G lies entirely in the original (flat) qq plane, then the spanning surfaces remain flat and the sum of areas 
A{1G) + A{G2) is RT, so the potential is zero. But for small transverse excursions of G into the x, y directions, the 
sum of areas is slightly increased. The simplest formula for this small change in the minimal surface area involves 
presenting the coordinates perpendicular to the surface, which we take to be x, y, as harmonic functions of two surface 
variables, which we take to be z, t: 



AilG) 



A{G2) - RT ~- I dtdz[ 

2 Jt(1G)+T{G2) 



(17) 



where the dot indicates a t derivative, the prime indicates a z derivative, and the integral is over the flat surfaces 
r(lG), r(G2) spanned by the original quark lines and the projection of the gluon contour G onto x,y — 0. At 
first sight this is confusing, because the gluonic path variables are originally presented as a Euclidean four-vector 
z^(t) = {x,y,z,t), with all four coordinates depending on the proper-time varible r. Because the gluon is assumed 
non-relativistic, Euclidean time t is essentially the same as proper time r. And the functions x(z,t),y{z,t) required 
in the area integral of Eq. H17|) are harmonic and therefore, in principle at least, determined by their boundary value 
on the contour (r). 
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Note that the two transverse variables x, y are just the variables used in to describe fluctuations of a surface 
formed from a chromoelectric flux tube extended in time. There should be, then, a close relation between the dynamics 
of a valence gluon in a hybrid and surface fluctuations. We will not explore that relationship further here, but it is a 
very interesting subject. 

Next, we assume that the actual gluon paths z^j^[t) stay close to the original minimal surface and arc not highly- 
convoluted. In that case, the perturbed surface is nearly flat. An example is a portion of the hclicoid, a minimal 
surface generated by rotating a rigid straight line segment around a fixed axis while at the same time moving it along 
the axis. Let the axis be the segment of the quark line 1 in the t direction, and the line segment at t = is the z axis 
from the origin to line 1. This line segment has length R/2. The line segment when t has a value scaling with the 
overall time T is slightly rotated around line 1 so that it passes through the point a: = 0, z = 0, y with y ^ R/2. The 
slightly-curved contour z^{t) is implicit in this description, and we need not be explicit; it is enough to observe that 
only small errors are made by approximating all portions of the contours by straight lines. Because the contour goes 
from a;, 2/ = at t = to small values y at t ~ T, it is clear that x ~ x/T, y ~ y/T. Since T ^ cxd we can drop the 
time derivatives. By a similar line of argument, x' ^ x/ R, y' ^ y/R- Then the area integral in Ea. p7|l is of the form 

A{lG) + A{G2)-RT^T[^{x^ + y^)]. (18) 

We need not work out the full details, because it should be clear by now that if we approximate pieces of the full 
contour by straight lines we are simply describing physically the gluon as a bead that runs without friction along 
a string stretched from z = R/2 to z = —R/2. Transverse string fluctuations stretch the string and give rise to 
the potential of Eq. H19|l below. As is evident from Eq. (|18ll . the small- fluctuation potential is that of a harmonic 
oscillator. 

Before discussing this potential (and its qqqG analog) we return briefly to the perimeter terms that keep the gluon 
from going too far beyond the original contour 12. These perimeter terms in the adjoint Wilson loop are discussed in 
the center- vortex picture in |14|.|15|. These terms give rise to a breakable string: A potential which is approximately 
linearly-rising, but goes flat when enough energy is stored in the string to pop a gluon pair out of the vacuum. It 
has been argued (see fl^ and references therein) that Ka is related to Kp by Casimir scaling; another elementary 
argument is that Ka = 2Kp because an adjoint string is two fundamental strings. The rough numerical estimates of 
[l5| yield a Ka that is bigger than Kp, but not quite according to either of these arguments. In the present work it 
is good enough, within the context of our approximations, to take Ka = 2Kp. 

Note that the adjoint string necessarily breaks when an energy comparable to 2M is stored in it. This means that 
the more relativistic the gluon, the closer one approaches to string breaking and mixing of the hybrid with states 
containing another gluon. We will not take this issue up here. 

III. GLUONIC POTENTIALS 

Rather than continue to deal directly with fluctuations of a minimal surface, as in Eq. p8|l . we will flnd a very good 
approximation to the potentials we need by starting from linearly-rising potentials attaching the gluon to the quark 
lines. Of course, these potentials cannot be taken literally when the gluon fluctuates far from the original minimal 
surface, and only their expansions for small fluctuations are of interest. As we saw in the last section, these give rise 
to harmonic oscillators with frequency ~ [Kp/{MR)Y/'^. Small fluctuations mean small gluon velocities, and so this 
regime of small fluctuations should be non-relativistic, if the gluon mass is large enough. 

Pot the qqG hybrid the gluon can move freely along the string, except that it cannot pass beyond the string 
endpoints. The free longitudinal momentum is associated with a momentum p in this surface of the form p — Nt:/R. 
If A'^ > 1 the motion is relativistic for R is less than about 2 fm, and so we will concern ourselves only with the 
ground state A'^ = 1. Details of an approximate non-relativistic analysis of the combined longitudinal and transverse 
Schrodinger equation are in the Appendix. 

As we will see, the numbers turn out such that even for the gluonic ground state of a hybrid relativistic corrections 
seem to be important, and we must ask where and how relativistic effects enter. We find in the Appendix that the 
gluon velocity P obeys ~ [A'i?/(M'^i?)]^/^, and relataivistic effects are important for R ~ 0.6 fm. 

Although we do not study mesonic hybrids in any detail in this paper, it is worth giving the general scheme of 
the gluonic potential for the qqG system, both for future use and because the same themes come up for the qqqG 
potential. 



FIG. 5: The potential of Eq. CHJi with R and Kf scaled to unity. The y=0 line is flat from a; = —1/2 to a; = 1/2. 

A. Mesonic hybrid potential 

For the qq system the arguments of the previous section lead to a non-relativistic potential that is a standard 
linearly-rising potential for each of the areas IG and G2 in Fig. ^ 

V{R, r)^KF{[(z- i?/2)2 + p2] 1/2 + [(^ + ^/2)2 + p^] 1/2 _ R] (19) 

[This potential has been proposed before .16] on purely phenomenological grounds.] This potential has all the prop- 
erties expected from our discussion in the previous section. In the minimal surface of the qq system, which is 
/9 = 0, \z\ < R/2, the potential is identically zero. For \z\ > R/2 there are linearly-rising terms beginning at \z\ — R/2 
and extending outside the minimal surface, as explained in the Appendix. We identify these terms with the break- 
able adjoint string, although the potential of Eq. I|19|) does not explicitly show breaking. For small excursions in p 
(transverse to the minimal surface) we expand the potential in p. The linear terms vanish, as they must (except at 
the endpoints \z\ = R/2), and near 2; ~ the quadratic terms are of the form (2Kpp'^ / R, with an oscillator frequency 
n(R) ~ 2[KF/iMR)]^/^ , where R is the qq separation. As z approaches the endpoints ±i?/2, the energy needed for 
a given transverse displacement grows, so there is a force keeping the gluon "bead" away from the endpoints unless 
the string on which the bead rides is not stretched at all. In Fig. [Slwe plot the potential for the d — 3 case, with q 
and q separated by unit distance along the x axis and y as the transverse coordinate. 
For further discussion of this potential, see the Appendix. 

B. Baryonic hybrid potential 

In general, the potential is the string tension Kp times the length of some minimal string network joining the original 
three quarks to the qq pair of the gluon. Wc will sec that in principle this potential requires consideration of minimal 
string lengths with different functional forms in different regimes. Once the work of expanding the potential in various 
regimes is all done, we get a result which differs only slightly from the result gotten from the intuitively-appealing 
potential 

V{Ri,R2, R3; r) = Kf{\Ri -r\ + \R2-r^ + {Rs - r^}- (20) 

The reader who is willing to believe this can go straight to the end of this section, where the final results for the 
harmonic terms are given. We will only work out the case when the three quarks (at fixed tirne) are at the corners 
of an equilateral triangle. This makes little difference in comparing theory with data, since |a] the corresponding 
lattice data are not much different for different quark triangles, provided that they have the same length L„iin of the 
Y-shaped Steiner string network joining them. 
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1 1 

(a) (b) 

FIG. 6: A view of the qqqG system looking along the quark world lines. The original Steiner network for no gluon is shown as 
dotted lines meeting at 0. (a) The gluon lies along the line 03; there is no resulting potential (except at the end points), (b) 
The gluon lies in regime 1 (defined in the text); the original Y strings are deformed into a longer Steiner network, raising the 
potential from zero. 

Fig. 1^ shows a view of a qqqG system looking down along the t axis on the transverse plane, taken to be the xy 
plane. The original quark system 123 has no angle as great as 27r/3 and so has an internal Steiner point; this is the 
only type of configuration we will consider here. Full circles represent the quarks 1, 2, 3 of Fig. |5]and the quark in 
the split gluon line, and the unfilled circle is the antiquark of that line. Dotted lines represent the original Y-shaped 
strings joining the quarks in the absence of the gluon; generically they meet at the point in Fig. Efa) at an angle of 
27r/3. Splitting the gluon leads to a new baryon 12G, which is joined by its own string network. 

We have chosen to single out quark 3 for special treatment in forming the baryonic hybrid. This could be appropriate 
if the quarks had different flavors or masses. Later, we will average the potential over interchanges of the three quarks 
to yield an energy symmetric in the quark labels. This is equivalent to dropping certain I = 2 terms in the potential. 

The essential point to keep in mind for the following discussion is that if a triangle has no angle as great as 27r/3, 
it has a Steiner point inside the triangle where the three strings meet at an angle of 27r/3. But if an angle of the 
triangle is 27r/3 or greater, the minimal-length string lies along the two sides meeting at that angle. Because of this, 
the algebraic form of the potential changes, depending on which of several regimes the gluon occupies. Fluctuations 
of the gluon in the out-of-plane (z) direction decouple from fluctuations in the XY plane, and are easily treated 
separately. If the gluon is fluctuating in-plane, we need to consider two regimes. Regime 1 is the sum of triangles 103 
and 203 in the figure; regime 2 is the triangle 012. 

Consider first regime 1. For the special position shown in Fig. EJa), with the gluon sitting on the original Steiner 
line 03, the strings for the new baryon coincide with those of the original. This is because the angle 102 is precisely 
27r/3. This, plus our remarks already made for a gluon lying in a mesonic surface, show that there is no potential 
energy associated with moving the gluon along the line 03, except within a QCD length of the ends and 3. Let us 
take the line 03 to lie in the y direction; then we expect to (and do) find only oscillations in the x direction. 

The potential to be expanded for all of regime 1 (again referenced to the original qqq potential) is 

V{Ri, R2, Rs'; r) — Kp{L„iin{RiG, R2G, R12) + I-R3 ^ L,nin{Ri2, R13, R23)} (21) 

where Ri is the vector position of the i*'' quark, p is the in-plane position of the gluon, Rij is the distance between 
quarks i and j, and Ric the distance between quark i and the gluon. Also, Ljnm{ci,b,c) is the length of the Steiner 
network joining three points that form a triangle of sides a, 6, c. The expression for this minimal length is: 

Lrmnia, 6, c) = -^{{a' +b' + c2) + V3[4a%^ - {a' + - c^)^] V2}i/2_ (33) 
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We take p = (x, 0) and find the quadratic term in the potential 



(23) 



where Lmin is the minimal length of the Steiner network joining the original three quarks. [If we had expanded the 
intuitive potential of Eq. (|2(J|I we would have found an equal contribution for y oscillations.] We will replace this by 
an average over choosing any of the quarks as having the gluon antiquark joined to it, instead of quark 3 as in the 

figure. This just amounts to isotropizing the result of Eq. (|23|l in the x y plane, with a resulting contribution 

from regime 1 to the potential of 

Vi^l^[x^ + y'). (24) 

This isotropization drops some 1 = 2, m = 2 terms in the potential. 

Regime 2 is the triangle 012 of Fig. El^a). If the gluon fluctuates (in-plane) into this triangle, the baryon triangle 
1G2 has an angle greater than 2tt/3. As a result, the minimal network is the sum of strings running from the original 
quarks to the gluon, and the potential is indeed that of the simple intuitive form in Eq. I|20() . This potential has the 
property that it is symmetric with respect to interchange of quark labels. It has an equilibrium point at the origin, 
reflecting the nature of the origin as a Steiner point. 

We give the expansion of the regime-2 potential for generic interquark distances, for future studies. Choose the 
(two-dimensional) quark vectors as 

Ri = ^2 - i?2(^, ^); i?3 = RsiO, 1). (25) 

With this coordinatization the Steiner point is at the origin. The expansion of the potential of Eq. 1)21(1 in transverse 
coordinates is 

l/.(i?i, i?„ i?3; P) = ^{p' E i" - E + • ■ • (26) 

where p = {x,y). Working out the algebra, using Eq. H25() . yields 

V2{RuR2,Rz]r) = V^, p, (27) 

where the matrix Vij has eigenvalues A^: 

\u \2 = -{y — ^-Y -^Y'^}. (28) 

This energy e(i?i, R2, R3) is clearly not a function of the minimal distance Lmin = Ri+ R2+ Rs- The deviations from 
dependence only on Lmin are acceptably small as long as none of the Ri are too small, which in any case is necessary 
for all of our developments. For the equilateral triangle, the regime-2 potential is 

V2 = l^{x' + y'). (29) 



Finally, the out-of- plane oscillations come from expanding Eq. H2U|) for a gluonic excursion in z. The result is 

V. = lJ^z^. (30) 



To construct the final potential, we note that regime 1 is twice the size of regime 2. For a simple harmonic-oscillator 
wave function as a variational trial wave function this suggests that we take the in-plane potential as 

V.n = lv, + \v2 = l^{x'+y'). (31) 



The total potential is 

V,n + V, = V{x) = '-^z'+'-^{x'+y') (32) 



PzKf 2 , PpKp , 2 , 2\ 
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FIG. 7: Lattice determination of the single-gluon excitation energy AE of a heavy qqqG system, as a function of Lmin, the 
length of the Steiner path joining the three quarks. From _5] . 



where (3^ — 9/2, (3p = 3/2. This potential can be isotropized in z but it makes little numerical difference. 

Strictly speaking, the harmonic potentials in V terminate at the positions of the quarks in the original baryon, 
and are replaced by something more complicated (essentially the breakable adjoint string). However, in the regime of 
large interquark separations this is not important. The reason is that if the unadorned harmonic potentials are used, 
the (non-relativistic) range of the gluon wave function is of order [Lmin/{MKF)]^^^- This is small compared to Lmin 
when Lmin ^ a, with a — {2KpM)~^/^ ~ 0.33 f. Even if this termination of harmonic potentials is taken into effect, 
it is a small effect numerically. 

The non-relativistic ground-state energy of the total potential V, plus the gluon mass M, yields a string energy for 
the baryon hybrid of 



3 LC p 

£string{Lmin) = M + {— + VS) [ 



1/2 - ^r ^1/2 
2 ' '^"MLmrn M L 



^-3.23. 



To this we add the color-Coulomb energy difference, estimated in the Appendix as 

V3 6 



Vc=ac{- 



2 -I- Qa2 

min ^ 



} 



(33) 



(34) 



with ttc — 0.15 5] and a, defined above, is ~ 0.33 f. (The length a accounts for spreading of the gluon wave function, 
as described in the Appendix.) 



IV. COMPARISON OF THEORY WITH LATTICE DATA 

There are extensive heavy-quark lattice data available on the heavy-quark qqqG system ^ . Fig. |3 shows the 
measured gluonic excitation energy AE of the gluon in the qqqG system, measured relative to the qqq system, and 
displayed as a function of the length of strings joining the original three quarks j5j]. Although in principle the energy 
depends on three separate interquark distances, @ finds that their lattice results fall pretty much on a single curve, 
parametrized by the minimum length Lmin of the Y-string joining the original qqq system. Both the energy and 
distance axes are scaled in lattice units, which for is about 0.1 fm. The energies in GeV are approximately twice 
the numerical values shown on Fig. [7|and on the next figure. 

We can match the data of Fig. \7\ reasonably well, using the variational treatment of the anisotropic harmonic- 
oscillator Schrodinger equation with relativistic kinetic energy, given in the Appendix, and a gluon mass M of around 
500 MeV. Coulomb effects are added perturbatively, as described in the Appendix. Since the data themselves fall 
roughly on a single curve for various qqq geometries with the same Lmin, we have done numerical calculations only 
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FIG. 8: Data (again in lattice units) of Fig. Qshown with the theoretical calculation discussed in the text. 

for the special case of an equilateral triangle. These calculations and the lattice data are shown in Fig. Ofor the 
value M = 500 MeV. 

The fit is surprisingly good at > 4, which is a distance short enough for relativistic effects to be important. In 

fact, the total energy at this distance is about 3Af , roughly apportioned equally among mass energy, potential energy, 
and kinetic energy. For smaller Lmin there really should be a fair amount of mixing with states containing more than 
one gluon, such as a baryon plus a glueball. We do not know to what extent such mixing is or is not accounted for 
in the lattice computations. The mass value of 500 MeV is a little small compared to a typical lattice value jsf of 600 
MeV, but certainly falls within an acceptable range of determinations of A/. A proper handling of mixing of higher 
angular momentum states might raise the fitted gluon mass somewhat. 

V. CONCLUSIONS 

We showed that there is a straightforward path from the fundamental elements of center-vortex QCD to a description 
of the gluonic ground state of heavy-quark qqqG hybrids (and, for that matter, of qq hybrids, although we have not 
attempted any comparison with lattice data for these). In graphic terms, the result is that the valence gluon can be 
thought of as a massive bead more or less freely running along a string joining the quarks, and subject to transverse 
restoring forces as the string is stretched by gluonic fluctuations. 

The hybrid problem simplifies at large interquark separation R, which is the regime we treated explicitly. Although 
our treatment, based on approximations to a Schrodinger equation with fuUy-relativistic kinetic energy is not exact, 
it is not subject to the usual burden of purely phenomenological assumptions made about the nature of confining 
potentials and so forth. The fundamental elements used were: 1) Splitting of a gluon proper-time adjoint Wilson 
line into a pair of co-located and oppositely-directed fundamental Wilson lines, which join to the three other Wilson 
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lines in a specific way. 2) Recognition that for a massive gluon the equal-time Schrodinger equation applies for 
sufficiently large R, but that relativistic corrections to the kinetic energy are important for physically-interesting 
quark separations. 

We used the gluon mass M as a single parameter for fitting the lattice data of Ref. 5] , and found a pretty good fit 
for a mass of about 500 MeV. This mass is a little less than the usual value of about 600 MeV found on the lattice, 
but at this stage of the investigation seems quite satisfactory. 

One could, in principle, think of extending this work to excited gluonic states in hybrids. However, at any reasonable 
interquark separation, the excited states are relativistic, and there will be some mixing with states with more than one 
gluon. As discussed in the Appendix for qqG hybrids, if this mixing is ignored there should be several contributions 
to the gluonic energy, including a stringlike energy of the form Nn/R, as invoked by Il2l| for such hybrids, plus 
terms of order [Kp / R)^/^ and Coulomb terms. 

We noted that the treatment of the valence gluon by adjoint-line splitting suggested a possible path toward describing 
the fluctuations of surfaces spanning fundamental Wilson loops. Originally, these surfaces were envisaged (l^] as having 
an essentially tangible existence as chromoelectric flux tubes. But fluctuations of spanning surfaces in the area law 
is more difficult to think about in the center vortex picture, where the description of confinement begins with the 
topological statement that the area law arises from the fluctuations in phase factors (elements of the center group) 
associated with topological linking of center vortices with the Wilson loop. At this level of description of confinement, 
the area in an area law is a minimal surface spanning the Wilson loop and is not subject to fluctuations However, 
coupling of gluons that are not part of the condensate of confining vortices can give rise to surface fluctuations, as we 
have seen in the present paper for valence gluons in hybrids. Work is underway to understand the description of such 
fluctuations when a large number of such gluons are coupled to a Wilson loop. 

Finally, we have not discussed here the mixing of the valence gluon in a hybrid with genuine qq pairs (as opposed 
to the fictitious pairs produced by splitting an adjoint Wilson line). If the valence gluon in qqqG is replaced by a 
valence qq pair, the result is a special kind of pentaquark, or qqqqq system. This pentaquark is definitely not the 
same as the famous pentaquark [l3| , which is of the form uudds, where the antiquark cannot pair with any quark to 
form a flavor-neutral state such as a gluon. Nonetheless, the problem of baryonic hybrid mixing with pentaquarks is 
of some interest, and worthy of future study. 



Although the gluon has a mass M, it is not heavy enough to insure that the non-relativistic Schrodinger equation 
applies throughout the range of quark separations for which lattice data are availaable. We give first a simple 
variational treatment of the ground-state excitation energy associated with the harmonic potentials encountered for 
hybrids, including relativistic effects for the kinetic energy. Then we go on to sketch a treatment of the gluon as a bead 
sliding freely along a string, and confined harmonically in the transverse direction. This is the problem encountered for 
the qqG hybrid and for one of the strings in the qqqG hybrid. Finally, we briefly sketch a treatment for color-Coulomb 
corrections. 



where Pz = 9/2, (3p — 3/2 and Lmin is the minimum length of the Steiner network. The full Hamiltonian, including 
relativistic kinetic energy, is 
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APPENDIX A: TREATMENT OF THE SCHRODINGER EQUATION 



1. Relativistic harmonic oscillator for the baryonic hybrid 



The string potential for the qqqG state is the sum of Eqs. (|30l I31|) . and has the form 




(Al) 



H = ^/f + M^ + V{x). 



(A2) 



We take a simple normalized variational wavefunction ip: 



IT IT Z 



(A3) 
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and calculate the expectation value {H) = e{az, ctp): 

{H) = —^ Hdp^p^ rdp\^{pl+pl + M'rf'exp{~^~Pi} + ^{^ + ^). (A4) 



The variational parameters Op^^ are to be determined by minimizing e, which can only be done numerically. Of 
course, the non-relativistic limit is easy to do, yielding 

(H) =M + (-^^) + 2^]. (A5) 

In the relativistic regime, one finds 

e=|(|^)^/^ (A6) 

where /3rei is a weighted average of /S^ and /3p, approximately given by /3rei = (l/3)(/3z + 2/3p). The separation between 
the relativistic and non-relativistic regimes is at a critical interquark distance i?c, where 

Rc-^. (A7) 

When R:^ Rc the problem is non-relativistic. Typically Rc is perhaps 0.8 f. 

The fit to data of Fig. |H1 comes from a combination of numerical and analytical work on minimizing expressions 
such as in Eq. ljA4|l . It turns out to be accurate enough to isotropize the potential of Eq. (|A1I setting a;^, and 
all to ?'^/3, and also to take az — ctp = a. Then a good analytic fit can be made to the kinetic energy term, and the 
result is 



3I3Kf 

t:ya) 's:^ [70; -|- ivi j ' + 

with 



£(«)^[7« + Af^]i/^ + ^ (A8) 



7~ 1.38, /3~ 2.52. (A9) 

The values of 7 and (3 are correlated so as to give the exact non-rclativistic energy at large R. It is now straightforward 
if lengthy to find an analytic expression for the minimum over a of e{a), involving the root of a quartic equation. 

2. Line of equilibrium 

We continue with a treatment of the Schrodinger equation that applies when there is a line of equilibrium, not just 
a point. 

To be specific, consider a a qq state, where the potential has the form in Eq. (|19|) (repeated for convenience) 

V{R,r) = Kf{[{z - R/2f + p^] V' + [(^ + R/^f + pY''' - R} (AlO) 

. Here (see Fig. 1) i? = (0, 0, R) is the separation vector of the heavy quark and antiquark, r — (x, y, z) is the gluon 
position vector, and = + y^ is the (square of the) gluon excursion out of the minimal surface. As discussed 
earlier, this is not literally the true potential, but it is accurate enough for our purposes. 
There is a second scale length a in the problem, given by 

a = {2KfM)-^/^ ~ 0.33/to, (All) 

and we will require i? 3> a. Since Rc/a ^ Kp/M"^ < 1, this requirement puts us in the non-relativistic regime, and 
means that our considerations may hold for the ground state, but not necessarily for excited states. 

We will see that the scale length for gluon excursions away from the minimal qq surface is at most of order 
{Ra^y^^ <C R- Therefore we provisionally take p^/R^ as an expansion parameter for the potential. Such an expansion 
requires \z\ <^ R/2, which we assume for now. The first two terms of the expansion are 

^^^_K^^SK^^^ (A12) 
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There are no terms in the expansion of the form 2:^; all powers of z are accompanied by powers of p^. Note that if 
p ^ 0, the point z = is the equilibrium point, midway between the quark and antiquark. 

The expansion above fails for \z\ > R/2. To characterize this regime we set p = in the potential, and find that the 
potential identically vanishes for \z\ < R/2, while for \z\ > R/2 there are linearly-rising potentials which confine the 
z motion. We will adopt the following essentially variational strategy, accurate enough for present purposes: Assume 
the ground-state wave function ip{r) for the Schrodinger equation is separable, of the form 

^^0(^)(^)i/2exp[-^] (A13) 

where ^(-z) is a normalized solution to the Schrodinger equation at p = 0. The other factor is a harmonic-oscillator 
ground state, and we will determine the parameter a from the potential of Eq. (jA10|l . averaged over the z motion. 
That is, the effective potential for p is given by 

y^ffip) = dz|0(z)|V(p,z). (AM) 

J-R/2 

We now construct a good approximation to (t>{z). For z > the potential is accurately given by 

V{z) ~ 2Kf\z ~ {R/2)\0{z - (i?/2)), (A15) 

and the region z < contributes symmetrically. The coefficient 2Kp is a decent approximation to the true adjoint 
string tension. Of course, the adjoint string is breakable, and the potential of Eq. I|A15|I is applicable only for 
sufficiently small |z — (i?/2)|. The ground-state solution to the Schrodinger equation which is symmetric under 

z — z is 

(/)(z) = Ni cos{qz)e[{R/2) - \z\] + {N2Ai{^)9[z - {R/2)] + z ^ -z}. (A16) 
Here A^i_2 are normalization constants, q is the momentum, and Ai(^) is an Airy function of argument 

e = i(z-f)-<zV. (A17) 
a 2 

The z-motion energy eigenvalue is Ez — /{2M), with the momentum q determined by the matching condition 

,qR Ai'iO ^ 
2 > 

At i? = this yields the well-known result (ga)^ ~ ^0 where —^0 — —1-02+ is the first zero of Ai'{£/). For large R one 
easily sees that q must decrease toward 7r/i?, so that Ez — tt^ /{2MR^). Analysis which we do not give here suggests 
that the following form should be a good representation of the momentum q{R) and Ez{R) as functions of R: 

'^^^^ - S^R-^ - 2Mam^l,R/af ^^^^^ 

1 /2 

where 7 = /n. This correctly yields the large- i? behavior and is correct to about 20% for small R. 

After some uninteresting algebra, one finds that at large R the Airy functions' renormalization constant N2 is small 
of order a/R, and therefore in calculating the p potential from Eq. I|A14|I we need the integral 

Veffip) = Kf dz^ cos2(gz){[(z - R/2f + p^Y'^ + [(z + R/2f + p^Y'^ - R}. (A20) 

J-R./2 

We have calculated this integral numerically, and find that for small {p/ R)'^ it is reasonably well-approximated by 
the z = expansion already given in Eq. HA12|I : V^ff — 2Kpp^ /R. The reason is that in the integral the weight 
function cos^ {qz) is small near the end points. Now the parameter a of the transverse wave function in Eq. ljA13|) is 
determined (a — 2{KpM/ Ry/^), and the total non-relativistic energy at large R is 

For this oscillator one has (p^) ~ (Ra^Y^^ which, as mentioned above, justifies expanding the full potential in powers 
of (p/i?)2 at large R. 

Excited-state solutions to the gluon motion in the z direction have momentum q given by ttN/R. When N > MR/tt 
the z motion is relativistic, with energy approximately equal to q. The transverse motion may be relativistic as well, 
requiring consideration of string breaking by gluon-pair formation. Without taking this phenomenon into account, 
one may expect corrections to the string energy Ntt/R to scale [as given in Eq. ljA6p ] with R~^/^. However, many 
other effects come into play, and we have not attempted to analyze data, such as that of (12j| . 



gatan(V) - -r-77rl?=-(9a)^- (^18) 
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3. Color-Coulomb corrections 

We will treat the Coulomb potential as a perturbation to the Schrodinger equation. If this potential were uniformly 
of order ad R for typical interquark separation i?, treating it as a perturbation would be valid (non-relativistically) 
when 

R » (A22) 

Kf 

which is easily satisfied for R values of concern to us. 

a. Baryonic hybrid 

In the baryonic hybrid qqqG, the minimum of the gluon potential is at the Steiner point, so it does not get close to 
the quarks, at least when these are well-separated (i? ^ a). When the quarks get close, so that R < a, the potential 
ceases to be 1/R, due to smearing by the gluon wave function. The net color-Coulomb potential for the quarks, which 
can be treated classically, in the original baryon is to be subtracted from the Coulomb potential for the hybrid. To a 
good enough approximation, the net Coulomb potential for our standard equilateral triangle is 

K = «c{7^-^=^=^=— }• (A23) 
We will use ac — 0.15, appropriate for the lattice with /3 = 6.0 and quenched quarks. 

b. Mesomc hybrid 

In the mesonic hybrid the gluon can get close to a quark even when these are well-separated, and the above argument 
does not straightforwardly apply when the gluon position z aliong the string is near ±i?/2, where the distance scale 
in the gluon-quark Coulomb potential is not generically of order R. Nonetheless, we will find that the probability 
that the gluon approaches the ends of its string, where the Coulomb potential is not of order ac/R, scales with 1/R 
so that in the end the Coulomb contribution to the energy does scale with 1/R. 

The Coulomb potential Vc (relative to the qq potential with no gluon) is 

Vc = " [(^ " ^/2)' + pT'^' - [(^ + Rm' + pT'^'}- (A24) 

A reasonable approximation is to replace the terms zzt (i?/2))^ in Vc by an expectation value of order a^, to replace 
by an expectation value of order [R/ {MKp)Y^'^ ~ a^^R/a)^^"^, and then to multiply the result by the probability 
that the gluon is within a distance a of z = ±i?/2. This probability is of order a/R at large R, as one sees from the 
matching of the wave functions given in Eq. (|A18|1 . Of course, when z <^ R the Coulomb potential is of order ac/R- 
The upshot is that all terms in the Coulomb potential are of this order: 

y -Aiac A2ac Oc > 

^~ i?[l+(il'/a)l/2]l/2 [^2+„2(^/„)l/2]l/2 + ^- ^ > 



where Ai_2 are positive coefficients of order unity. Since we are not treating the mesonic hybrid in detail, we have 
neglected to write some other factors of order one at various other places in this equation. 
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